Freezing of spin dynamics in underdoped cuprates 
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The Mori's memory function approach to spin dynamics in doped antiferromagnetic insulator combined with 
the assumption of temperature independent static spin correlations and constant collective mode damping leads 
to uj/T scaling in a broad range. The theory involving a nonuniversal scaling parameter is used to analyze 
recent inelastic neutron scattering results for underdoped cuprates. Adopting modified damping function also the 
emerging central peak in low-doped cuprates at low temperatures can be explained within the same framework. 
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It is experimentally well established that magnetic static 
and dynamical properties of high-T c cuprates are quite 
anomalous. Even at intermediate doping the magnetic re- 
sponse does not follow the Fermi-liquid phenomenology in 
the metallic state, and reveals the remarkable resonant peak 
within the superconducting (SC) state. With decreasing 
doping the magnetic response approaches in a novel way 
the one of the reference undoped material-the antiferromag- 
netic (AFM) insulator. It has been first observed by in- 
elastic neutron scattering (INS) experiments on low-doped 
La 2 - K Sr x Cu0 4 (LSCO) HI Hfl that local (q-integrated) spin 
dynamics in the normal state exhibits anomalous uo/T scaling, 
not reflected in the AFM correlation length £ which shows no 
significant T-dependence. Similar behavior has been found 
also in YBaCu 3 6+;l; (YBCO), initially at x = 0.5, 0.6 H 
with the regime restricted to T > T c , and in Zn-doped YBCO 
H with reduced T c ~ 0. 

Recent, more detailed INS experiments on heavily un- 
derdoped (UDX including Li-doped LSCO |H, YBCO with 
x = 0.35 0, 0] and x = 0.45 fl EJ, as well as electron- 
doped Pri_ x LaCe x Cu0 4 _,5 (PLCCO) with x = 0.12 ifloiflltl, 
confirm universal features of anomalous normal-state spin dy- 
namics. Latter materials exhibit no SC (Li-doped LSCO) or 
have very low T c (YBCO, PLCCO) so that the lu/T scaling is 
found in a broad range both in q-integrated local susceptibility 
Xl( w ) lUlSS and in ViUo) at the commensurate AFM q = 
Q = (tt, tt) USEE]]]. Typically, x'qM is a Lorentzian 
with the characteristic relaxation rate scaling as Y — aT, but 
with a nonuniversal a. Similarly, normalized x'l( w ) nas been 
fitted to the scaling form f(u,T) = x£(w, T)/x'l(u, T = 
0) ~ (2/7r)tan- 1 (Aiw/T + A 2 (u/T) 3 ) fiH with material 
dependent A\2. It is characteristic that the normal-state-like 
scaling is found even below T c without any clear sign of the 
presence of the resonant peak. On the other hand, it has been 
observed that at low T the spin dynamics gradually transfers 
into a quasi-elastic peak for T < T g (J, |6(] whereas the in- 
elastic response saturates. However, this freezing mechanism 
appears to be of entirely dynamical origin since £ as well as 
the integrated intensity are unaffected by the crossover. 

In order to describe the anomalous scaling of spin dynam- 
ics as well as of other electronic properties of cuprates the 
concept of marginal Fermi liquid has been introduced 111 211 . 



Whereas the deeper origin of such behavior is often given 
in terms of proximity to a QCP 11131 Il4tl the above systems 
clearly lack the criticality of £(T). It is thus quite remark- 
able that a similar behavior is also found in heavy-fermion 
metal system CeCuAu H15I1 which possesses a well defined 
QCP (unlike in cuprates) in the phase diagram, nevertheless 
exhibits similar uj/T scaling without the criticality in £(T). 

The present authors introduced a theory of spin dynam- 
ics in doped AFM IU6I1 which describes the scaling behavior 
as a dynamical phenomenon. Namely, assuming that static 
quantities-the equal-time spin correlations in particular-are 
unaffected (or weakly dependent) by T, and that the system is 
metallic with finite spin collective-mode damping, the system 
close to AFM naturally exhibits ui/T scaling in a wide en- 
ergy range, but unlike at usual QCP, with a saturation at low- 
enough T. It should be noted that such a scenario is close to 
the theory of freezing of a liquid, which is as well a dynamical 
phenomenon reflected only weakly in equal-time correlations. 

The aim of this paper is to analyse recent INS experiments 
on different cuprates mentioned above within this framework 
and show that observed differences (nonuniversality) in scal- 
ing behavior can be made consistent taking into account actual 
parameter regimes. Moreover, the analysis can be generalized 
to the emergence of the CP assuming same static quantities 
but the change in collective-mode damping. 

Within the memory-function approach |3 17] the dynam- 
ical spin susceptibilty can be generally expressed as 



UJ 2 + CjM q (w) — 



(1) 



The spin stiffness r/ q is the first freqency moment of the 



dynamical susceptibility r; q 



[Sta,S')]}, and can be 



within any relevant model expressed in terms of equal-time 
correlations only, whereas w q = (f/q/Xq) 1 ^ 2 is an effective 
collective mode frequency where x q = Xq( w = 0) is the static 
susceptibility.Since we have in mind a doped AFM the actual 
evaluation in relevant models like the t-J model JlH 17 1 re- 
veal that ?/ q is only weakly q-dependent for q ~ Q whereby 
?7q ~ '7 ~ 

The form (UJ is particularly suited for a description of 
damped collective modes close to an ordered magnetic state. 
Mq(cj) is the (complex) memory function incorporatinginfor- 
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mation on the collective-mode damping 7 q (w) = M''(to). 
Since we will be here dealing withdoped cuprates in the nor- 
mal state which are (anomalous) paramagnetic metals (PM), 
low-frequency collective modes at q ~ Q are generally over- 
damped 7 q > uiq consistent with INS results. 

Underdoped cuprates close to the AFM phase represent a 
system of low charge-carrier concentration but large spin fluc- 
tuations whose dynamics is quite generally restricted by the 
sum rule 



results in 



1 



cth— ^(w) = (Si 



(2) 



In doped AFM C q is strongly peaked at Q with a character- 
istic width kt = Moreover, the total sum rule is for 
a system with local magnetic moments (spin 1/2) given by 
{ 1 /N) y"L Cq = (1 — c/j)/4, where Ch is an effective hole 
doping UJ]. 

While the formalism so far is very general, we now intro- 
duce approximations specific to UD 11611 . INS experiments 
listed above indicate that within the normal state the effec- 
tive q width of Xq( w )> i- e -> dynamical k(u>), is only weakly 
T- and cj-dependent, even on entering the regime with the CP 
response which implies that C q is T independent or at least 
not critical for T — » 0. This has been also confirmed be the 
present authors [16] in numerical investigations of the rele- 
vant t-J model, where kt has been found to be rather large 
even for the lowest doping. Since ?? q ~ rj is also nearly 
constant the to/T scaling and the freezing mechanism must 
emerge from spin dynamics not reflected in equal-time corre- 
lations. Within further analysis we will assume the commen- 
surate AFM response at Q and the double-Lorentzian form 



<7 q = C/[(q-Q) 



„2 12 



close to INS experiments although 



qualitative results at low to do not depend on a particular form 
ofC q . 

Paramagnetic metal: Even at rather low doping cuprates be- 
have as PM, as manifest, e.g., by their metallic resistivity 
p(T) and quasiparticle-like excitations on (the parts of) the 
Fermi surface as revealed by ARPES. This is underlying our 
assumption that the (Landau) damping 7 q (oj) is dominated by 
particle-hole excitations [17] being weakly q and to depen- 
dent, hence 7 q (w) ~ 7, as has been found also in the nu- 
merical analysis of model results 11611 . whereby a noncritical 
dependence j(T) can still play some role. We treat here 7 as 
a phenomenological parameter although it can be estimated 
from model results 111611 or even analytically 11711 . It could 
be also extracted from high-energy INS experiments follow- 
ing the crossover to underdamped modes tOq ~ 7. How- 
ever, most INS data on UD cuprates are restricted to energies 
to < 50 meV being too low to observe the crossover. 
The assumption of constant 7 in Eq. ([]]) leads to 



TfftO 



{to 2 — to" 2 .) 2 + 7 2 w 2 



(3) 



L0T r 



7Xq' 



(4) 



While such form appears plausible, it can be explicitly tested 
experimentally. Recent INS data obtained for Xq( w ) have 
been fitted with Eq.© [5, 7] to extract Tq(T). Indeed, plot- 
ting the interdependence of data for Xq and Tq reported in 
10] (see Fig. 12) shows that Xqr Q (T) = const. (= 77/7) to 
within experimental error over the whole span of Tq, confirm- 
ing in this way our assumption of constant 7(0;) ~ 7. 

Let us now concentrate on the behavior of Xq(w) as en- 
forced by the sum-rule, Eq.©. It has been noticed that the 
variation of Tq is mainly determined by the parameter 



C = C Q 7r 7 /(2ry), 



(5) 



and only marginally on, e.g., 7 provided that we are discussing 
the regime T « 7 which is experimentally relevant. 




Figure 1: Relaxation rate Tq vs. T for different parameters £ and 
7=40 meV. Inset: dependence of slope parameter a on £ for two 
different 7. For comparison, an estimate for a based on Eq. {6} is 
included. 

In Fig. 1 we first present results for Tq(T) in the experi- 
mentally relevant regime and units for a range of £ = 1 — 8. It 
is evident that for any £ there exists a finite limiting value 
Tq = Tq(T — » 0) being strongly dependent on (, i.e., 
rgj - 7exp(-2C) [16]. On the other hand, for T > Tq 
the variation is nearly linear Tq ~ aT being a manifestation 
of the to/T scaling. It is also seen that the linear-in-T depen- 
dence of Tq extends well beyond the limit to ~ 7/2 provided, 
of course, that kt remains weakly dependent on T in the re- 
spective temperature region. 

As seen in Fig. 1, the slope a is not universal but depends 
on £ as well as on 7, where the latter dependence is very weak 
(logarithmic) as is also manifest from the collapse of the data 
for two values of 7 onto a single set (inset to Fig. 1). In fact, 



which in the overdamped regime, w q < 7/2 and to <C 7, from Eqs. I© and OJ and for w q < 7/2 one obtains 
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*(1 



J_Q_, 

2nT' 



-In 



7 1-1 



2ttT j 



(6) 



provided that 7 ^ 27rT, and ^f(-) is the digamma function. 
Hence, in the regime Tq ~ aT and for moderate a < 2ir, we 
get a ~ 7r/ (2£) within leading order, which is seen to agree 
very well with the exactly calculated values (inset to Fig. 1). 

Obtained results can be directly put in relation with INS ex- 
periments on UD cuprates. Several results have been recently 
reported for Xq(w), scanned in both q and u> for different T. 
In particular, it has been observed that Xq^) can be well de- 
scribed by the simple Lorentzian, Eq. and T(T) ~ aT has 
been extracted. For UD YBCO with x = 0.35 the authors [7] 
give a = 0.7, whereas for LSCO weakly doped with x = 0.06 
Li scaling is well described by a = 0.18 [5]. Electron-doped 
PLCCO with doping x = 0.12 also exhibits scaling ifTHl with 
a = 0.44. It is therefore plausible that different a can be 
within our theory explained with different £, Eq.©. Since 
Cq oc the strongest dependence results from kt and 

experimental results onctal/( are consistent with the rank- 
ing of measured (low-frequency) K, i.e. k ~ 0.03 (in r.l.u.) 
for YBCOJ0L_ k ~ 0.01 for LSCO H, and k ~ 0.02 for 
PLCCO H 1 OL 1 1 lfl . Since 77 is a robust quantity and rather well 
known from model calculations, only unknown 7 prevents so 
far more quantitative statements. 

The above experimentally deduced a all require quite large 
£, hence one can expect, following Fig. 1, very low satura- 
tion r q for T — > 0. On the other hand, mentioned INS data 
for YBCO and LSCO seem to indicate quite substantial Tq. 
However, saturation of Tq, setting in for T < T g , is accom- 
panied by the simultaneous appearance of the CP. This effect 
goes beyond here presented (simple) explanation and requires 
within the PM the introduction of additional damping, as dis- 
cussed lateron. 

Recent detailed INS measurements by Hinkov et al. on 
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Figure 2: Temperature evolution of x'q(u,T) for UD YBCO with 
x — 0.45 (symbols) compared with theoretical result, Eq. 10, 
with C, = 1.8 and 7=40 meV (lines). 



UD YBCO with x = 0.45 Mm allow also for a more com- 
plete analysis of Xq(w) as a function of T. Data normalized 
at T = 0, i.e., x'q(w,T)/x q (cj,T = 0) are presented in 
Fig. d2J and compared with theoretical ones, following from 
Eqs.ffSl),©, where the only relevant parameter is £. The over- 
all agreement is quite satisfactory and is achieved for ( ~ 1.8 
implying Tq ~ 1.3 meV. The marked dicrepancy between 
theory and experimental data at oj = 3 meV is tentatively as- 
sociated with a possible saturation of Tq (%q(w, T — 0)) to 
a larger (smaller) value then predicted by theory. Namely a 
better fit to experimental data based on £ and T° as indepen- 
dent parameters, yields unchanged ( but substantially larger 
Tq ~ 4 — 5meV, which would correspond to a crossover 
temerature T g ~ Tq ja ~ 50K. On the other hand, the poor 
agreement with the 32.5 meV data should be attributed to the 
breakdown of oj/T scaling since u> > 7/2. 
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Figure 3: Scaling of theoretical normalised x'1{uj,T) based on £ 
and 7 as in Fig. For comparison, the scaling function /(x) as 
frequently used in fits to experiments is also plotted. 

Experimental results for q-integrated Xz( w ) m the same 
regime also reveal the uj/T scaling and have been several 
times fitted with the ansatz involving tan _1 (a;), as already 
noted in the introduction. However, the fit did not find so 
far a deeper theoretical background. It has been nevertheless 
shown that results of the present theory are quite close to this 
description [16]. In fact, assuming a Lorentzian dependence 
of Xq on q and invoking the relation T q = r)/(jx%) a simple 
calculation yields tan _1 (Aiu;/T) with A\ ~ 1/a, provided 
that k 1 but k 2 Xq ~ const. Note, however, that quite gen- 
erally Xq obtained from Eq. (01 decays at least as a double- 
Lorentzian, so that Xl( w j T) is a more complicated function, 
allowing even for nonanalytic behavior in u> for £ ^> 1. 

In Fig. (3) we next present the theoretically obtained data 
for the scaling function /(co/T) = Xl( w > T = 0) 

for several different lu and the temperature range — 300 K 
based on the same set of data for 7 and £ as discussed in rela- 
tion to Fig. (2). While the co/T scaling is quite remarkable it 
breaks down at higher energies, as expected. We also plot the 
scaling function f(x) = 2/7rtan _1 (Aix + A2X 3 ) often em- 
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ployed in fits to experimental data: the 'best-fit' coefficients 
A\ ~ 0.56 and A^ ~ 0.09 are in reasonable agreement with 
values A\ ~ 0.49, A 2 ~ 0.12 determined directly from ex- 
perimental data ||9|] . 

CP response: The advantage of the memory-function for- 
malism is that the emergence of the CP at T < T g in the spin 
response can be as well treated within the same framework. 
One has to assume that unlike in a PM the mode damping 7 q 
is not constant but may acquire an additional low frequncy 
contribution. In particular we can take 



M a 



ij - 5 2 /{uj -HA), 



(7) 



with T-dependent S and A, which via Eq.© leads to Xq( w ) 
of the form used also to analyse experimental INS data for 
YBCO with x = 0.35 (g]. It has been first employed for the 
analysis of the CP in ferroelectric s [ 18], assuming a coupling 
of a damped mode to a slowly fluctuating object. 

For A — > but S 2 /\ 3> 7 the modified M q leads to two 
distinct energy scales and hence to two contributions to spin 
dynamics, i.e., the CP part x q ( w ) ar, d the regular contribution 



Xq( w ) 



x q r c 



Aq 1 r 

T r — U 



q ~ Q, (8) 



valid for u> < A and A < w « 7, respectively. Thus, below T g 
the new scales are set by T r — fi q /7 and r c = {j]/8 2 )\/x% 
with x q ° = v/^q an d ^ q = ^ q + S 2 . If one assumes 
that 6 saturates at low T, as is manifest by saturation of T r 
0], T c oc A/ Xq becomes the smallest energy scale, result- 
ing in a quasielastic peak which is, according to Eq. ([HJ, a 
Lorentzian of width T c , with the static value Xq( w ~~ v 0) = 
X q , as required by Eq. ([TJ. However, the real part drops 
for uj > A to the value x q °, hence the CP contribution is 

A Xq = Xq (W 2 «X q - 

It should be pointed out that within such a formalism a mod- 
ified Af q does not (necessarily) induce a change of equal-time 
correlations as C q , in particular the characteristic low-w width 
k, remains unchanged, a fact observed experimentally |^,|7|]. 
Moreover, as the CP is too narrow to affect the moments of 
Xq(u>), i]q is determined solely by the regular part. It is then 
plausible to assume that also x q "(w) remains constant on low- 
ering T < T g , as also observed experimentally in Li-doped 
LSCO [5] where below T g ~ 50 K the q-integrated regular 
part takes the form \'i ( w ) kw/ (lu 2 +T 2 ) with T-independent 
r. The CP contribution then follows from the sum rule 



TA Xq = C^T) = C q - C*(T), 
with the regular part given by [ IS 



CTJT) = -Mn 



77 7 ttT 



T<T„ 



(9) 



(10) 



7T7 r r jT r 

again consistent with experiment on x — 0.35 YBCO [7] 
where below T ~ 50 K the (integrated) intensity transfers 



from a broad inealstic scattering to a CP while conserving the 
(total) sum rule. 

From the above considerations it follows that the onset of 
CP can be explained as a dynamical freezing phenomenon. 
However, it is evident that the appearance of the anomalous 
(singular at low u>) damping M q (o;) is not consistent with nor- 
mal metallic particle-hole excitations. We can speculate that 
to obtain such a form one has to invoke an inhomogeneous 
(disordered) system with predominantly localized charge car- 
riers preventing low-w scattering, resulting in a nearly di- 
verging M q (w —►()). In any case the transition to CP-type 
dynamics can be associated with a significant softening of 
A(T < T g ). 

In conclusion, we have shown that the approach presented 
gives a consistent explanation of the uj/T scaling both in 
X q (w) as well as in x'[ It i s based on two well established 
experimental facts: the overdamped nature of the response 
and the saturation of effective inverse correlation length k(lu) 
at low u> and T leading to the anomalous slowing down, i.e., 
gradual freezing of the dynamic spin response, a phenomenon 
of entirely quantum origin. The theory can as well account for 
the INS results in UD cuprates, yielding an explanation for the 
nonuniversal parameter a = Tq/T. 

The appearance of the CP for T < T g is easily incorpo- 
rated into the formalism via the (almost) singular ui- and In- 
dependent damping M q (aj). However, the question as to the 
nature of CP, whether of intrinsic origin, based on spin-glass- 
like arguments, as suggested, e.g., in (01, or of some other, 
possibly even extrinsic origin, remains to be settled. 

The authors would like to thank Dr. V. Hinkov for using his 
data prior to publication and for fruitful discussion. This work 
was supported by Slovenian Research Agency (ARRS) under 
grant No. Pl-0044-106. 



[1] B. Keimer et al, Phys. Rev. Lett. 67, 1930 (1991); Phys. Rev. 
B 46, 14034 (1992). 

[2] for a review see M. A. Kastner, R. J. Birgeneau, G. Shirane, and 
Y. Endoh, Rev. Mod. Phys. 70, 897 (1998). 

[3] B. J. Sternlieb, G. Shirane, J. M. Tranquada, M. Sato, and S. 
Shamoto, Phys. Rev. B 47, 5320 (1993). 

[4] K. Kakurai, S. Shamoto, T. Kiyokura, M. Sato ,J. M. Tran- 
quada, and G. Shirane, Phys. Rev. B 48, 3485 (1993). 

[5] W. Bao, Y. Chen, Y. Qiu, and J. L. Sarrao, Phys. Rev. Lett. 91, 
127005 (2003). 

[6] C. Stock, W. J. L. Buyers, Z. Yamani, C. L. Broholm, J.-H. 

Chung, Z. Tun, R. Liang, D. Bonn, W. N. Hardy, and R. J. 

Birgeneau, Phys. Rev. B 73, 100504(R) (2006). 
[7] C. Stock, W. J. L. Buyers, Z. Yamani, Z. Tun, R. J. Birgeneau, 

R. Liang, D. Bonn, and W. N. Hardy, Phys. Rev. B 77, 104513 

(2008). 

[8] V. Hinkov, D. Haug, B. Fauque, P. Bourges, Y. Sidis, A. Ivanov, 
C. Bernhard, C. T. Lin, and B. Keimer, Science 319, 597 (2008). 

[9] V. Hinkov, Ph. D. thesis (2007), unpublished. 
[10] S. D. Wilson, S. Li, H. Woo, P. Dai, H. A. Mook, C. D. Frost, 
S. Komiya, and Y. Ando, Phys. Rev. Lett. 96, 157001 (2006). 



5 



[11] S. D. Wilson, S. Li, P. Dai, W. Bao, J.-H. Chung, H. J. Kang, 
S.-H. Lee, S. Komiya, Y. Ando, and Q. Si, Phys. Rev. B 74, 
144514 (2006). 

[12] C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. Abrahams, 

and A. E. Ruckenstein, Phys. Rev. Lett. 63, 1996 (1989). 
[13] S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. Rev. 

Lett. 60, 1057 (1988); Phys. Rev. B 39, 2344 (1989). 
[14] G. Aeppli, T. E. Mason, S. M. Hayden, H. A. Mook, and J. 

Kulda, Science, 278, 1432 (1997). 
[15] A. Schroder, G. Aeppli, E. Bucher, R. Ramazashvili, and P. 

Coleman, Phys. Rev. Lett. 80, 5623 (1998). 



[16] P. Prelovsek, I. Sega, and J. Bonca, Phys. Rev. Lett. 92, 027002 
(2004). 

[17] I. Sega, P. Prelovsek, and J. Bonca, Phys. Rev. B 68, 054524 
(2003). 

[18] B. I. Halperin, and C. M. Varma, Phys. Rev. B 14, 4030 (1976); 

S. M. Shapiro, J. D. Axe, G. Shirane, and T. Riste, Phys. Rev. 

B 6, 4332 (1972); I. Sega, Ph. D. thesis (1979), unpublished. 
[19] Note that this estimate follows from Eq. ([6} on assuming that 

r q (T) saturates to (almost) T-independent r r below T g . 



